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TITLE OF THE INVENTION 

DELAY TIME ESTIMATION METHOD AND 
RECORDING MEDIUM STORING ESTIMATION PROGRAM 



5 BACKGROU ND OF THE I NVE NT ION 

1. Field of the Invention 

The present invention generally relates 
to delay time estimation methods for estimating a 
delay time in a logic circuit composed of 

10 transistors and to recording mediums in which a 
computer program for executing a delay time 
estimation method is stored. More part icularly , the 
present invention relates to a delay time 
estimation method and a recording medium storing an 

15 estimation program with which efficient and precise 
estimation of delay time is possible. 

2. Description of the Related Art 

The necessity for accurate signal 
analysis in a logic circuit is growing for 

20 designers of high-speed large-scale LSIs. Signal 
delay is one of the most important parameters 
because it is important for an LSI designer to know 
a timing margin in order to determine whether an 
LSI can operate properly. Methods for modeling and 

25 estimating time delay have been proposed. 

Figs. 9 through 11 show a processing 
flow for delay time estimation according to the 
related art. Referring to Fig. 9, the related-art 
delay estimation includes extraction of information 

30 relating to connection of the target circuit from a 
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layout. The extracted circuit is modeled as a 
series comprising inverters (INV) 101 and 102 
connected via a wire 103. Based on this circuit 
connection information, circuit configuration 
5 information in which a load is modeled by an RC 
component is produced, as shown in Fig. 10. A 
source model 106 corresponding to the inverter 101 
is modeled by a combination of a power source 104 
and resistance 105. A load component 108 comprises 

10 an RC distributed constant circuit 107 

corresponding to the wire 103 , and input pin 
capacitance Cg corresponding to the inverter 102, 
the RC distributed constant circuit 107 and the 
capacitance Cg being determined so that the 

15 admittance downstream from an output terminal of 
the gate provides a match with a third-order 
approximation. 

The resistance and capacitance 
constituting the load component 108 is approximated 

20 by a finite number of RC components. The input 

capacitance of the RC distributed constant circuit 
107 is modeled by capacitance C2 . The combined 
output capacitance of the RC distributed constant 
circuit 107 and input pin capacitance Cg of the 

25 inverter 102 is modeled by capacitance Cl. The 

capacitance Cl, capacitance C2 and resistance R of 
the RC distributed constant circuit 107 form a jc- 
load model 109 as shown in Fig. 11. 

The jt-load model, formed by two C 
30 components and one R component, is constructed such 
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that, for any type of source model 106, a voltage 
waveform, occurring at the output terminal of the 
gate when the circuit of Fig. 10 is established, is 
approximated by a voltage waveform obtained as a 
5 result of analyzing the circuit of Fig. 11. 

The approximation described above is 
disclosed in Peter R. O'Brien and Thomas L. 
Savarino, Modeling the Driving-Point Characteristic 
of Resistive Interconnect for Accurate Delay 

10 Estimation, Proc. IEEE International Conference on 
Computer-Aided Design, 1989. Admittance Y(s) 
looking downstream from the gate output terminal is 
estimated (Fig. 3, equations (19)-(32)). The 
admittance Y(s) thus obtained is used to estimate R, 

15 CI and C2 (equations (14)-(16)). 

The load model 109 is then connected to 
the source model 106 so as to estimate a delay time 
through response analysis. The voltage level of the 
power source 104 and the resistance 105 have 

20 respective values determined by modeling conditions. 
The method of computing the voltage level and 
resistance is described in details in Florentin 
Dartu, Noel Menezes, Jessica Qian, and Lawrence T. 
Pillage, A Gate-Delay Model for High-Speed CMOS 

25 Circuits , Proc. 31 st ACM/ IEEE Design Automation 

Conference, 1994, so that a detailed description is 
omitted. 

A description will now be given of the 
operation according to the related art. 
30 Fig. 12 shows a construction of the 



inverter 101 comprising a PMOS transistor and an 
NMOS transistor. In a rising transition at the 
output terminal Y, as the potential at the input 
terminal A goes from a low level (L) to a high 
5 level (H), the PMOS transistor PI makes a 

transition from an OFF state to an ON state so as 
to charge an output load. When an increase in the 
potential between the source and drain of the PMOS 
transistor PI is relatively smaller than the 

10 magnitude of change in the gate potential, a 
transition from a region, characterized by an 
increase in a current with time, to another region 
characterized by a rapid exponential decrease in 
the current occurs (see pattern 2 of Fig. 5). 

15 Referring to Fig. 13, in the related-art source 

model 106, an internal voltage source E(t), whose 
voltage level shows a linear variation between 0 ad 
Vdd in a time At, is used to represent the 
transition described above. 

20 According to the related art, the ot-load 

model 109 shown in Fig. 11 subject to delay 
estimation is approximated by a purely capacitive 
load model as shown in Fig. 14 providing an 
equivalent response. In this case, equivalent 

25 capacitance is determined by considering the 
shielding effect provided by the resistance R 
constituting the jt-load model 109. A delay time is 
determined by searching a table listing delay time 
along with gradients of predetermined input 

30 waveforms and output load capacitance. The table is 
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searched so that a delay time that matches the 
modeling condition is determined by interpolation. 

The related-art delay time estimation 
has a disadvantage in that it is not adapted for 
5 another possible transition pattern (pattern 1 of 
Fig. 5) in which there is a transition from a first 
region, characterized by an increase in current 
with time, to a second region, characterized by a 
gradual decrease in current, and then to a third 

10 region, characterized by an exponential decrease. 

The related-art internal power source model E(t) as 
shown in Fig. 13, characterized by a linear 
variation between 0 and Vdd in a time At, fails to 
represent a saturation region in which the current 

15 gradually decreases (region 2 of Fig. 5). Therefore, 
the related-art method fails to provide delay time 
estimation that matches the operating 
characteristic of transistors. 

Another disadvantage of the related art 

20 is that there is a need for a library of two- 
dimensional delay tables listing gradients of input 
waveforms and output load capacitance, thus making 
it necessary to store a large volume of data. 
Interpolation errors are incurred as a result of 

25 using the tables. If the rc-load model is to be used 
instead of the purely capacitive model, the 
dimension of the table increases so that the volume 
of data is increased, thereby rendering its 
implementation impossible. The practice of 

30 conversion into equivalent capacitance, performed 
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in this background, generates errors. 

SUMMARY OF THE INVENTION 

Accordingly, a general object of the 
5 present invention is to provide a delay time 

estimation method and a recording medium storing an 
estimation program in which the aforementioned 
disadvantages are eliminated. 

Another and more specific object of the 

10 present invention is provide a precise delay time 

estimation method and a recording medium storing an 
estimation program in which a saturation region, 
characterized by a gradual decrease in current, is 
represented so that delay time estimation that 

15 matches the transistor characteristic. 

The aforementioned objects can be 
achieved by a delay time estimation method for 
estimating a delay time in a logic circuit that 
includes a MOS transistor, comprising the steps of: 

20 modeling the MOS transistor by a resistive element 
having fixed resistance and a power source voltage 
that varies with time; and segmenting an operating 
characteristic of the MOS transistor thus modeled 
into a fist region in which a current increases as 

25 a gate potential varies, a second region 

corresponding to a saturation region of the MOS 
transistor in which region the current gradually 
decreases as the gate potential remains constant, 
and a third region corresponding to a linearity 

30 region of the MOS transistor in which region the 
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current decreases as the gate potential remains 
constant. 

The delay time estimation method may be 
adapted for a circuit in which a plurality of logic 
5 circuits that includes MOS transistors and comprise 
the steps of: segmenting an operating 
characteristic of last-stage MOS transistor 
constituting a logic circuit of a last stage into a 
first region in which a current increases as a gate 

10 potential varies, a second region corresponding to 
a saturation region of the last-stage MOS 
transistor in which region the current gradually 
decreases as a gate potential remains constant and 
a third region corresponding to a linearity region 

15 of the last-stage MOS transistor in which region 

the current decreases as the gate potential remains 
constant. 

E=Rsxi(t)+v(t) may hold for t=At x and 
t=At x +At 2 , where E denotes the power source voltage, 

20 Rs denotes resistance of a model of the power 

source, i(t) denotes a charge current of a load 
model, v(t) denotes a charge voltage of the load 
model, and wherein V lf At x and At 2 are determined 
based on a fact that values of E-v(t) and i(t) 

25 reside on an Ids-Vds characteristic curve at a 

given gate potential, where Ids denotes a drain- 
source current and Vds denotes a drain-source 
voltage, and where V 1 denotes a voltage at a 
boundary between the first region and the second 

30 region, At x denotes a time required to arrive at 
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the boundary, and At 2 denotes time required to 
reach the power source voltage via the second 
region. 

The delay time estimation method may 
5 employ a delay library including function 

information for specifying polygonal lines that 
provide a model of an Ids-Vds characteristic at a 
given potential and also including function 
information related to a slew rate specifying a 
10 fixed delay. 

The aforementioned objects can also be 
achieved by a recording medium storing a computer 
program that executes a delay time computation 
method . 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other objects and further features of 
the present invention will be apparent from the 
following detailed description when read in 
20 conjunction with the accompanying drawings, in 
which: 

Fig. 1 illustrates extraction of circuit 
connection information according to a first 
embodiment of the present invention; 
25 Fig. 2 illustrates generation of circuit 

configuration information; 

Fig. 3 illustrates generation of load 

model; 

Fig. 4 is a graph showing variation of 
30 power source voltage with time; 
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Fig. 5 is a graph showing patterns of 
waveform of a current; 

Fig. 6 is a Ids-Vds characteristic graph 
according to the E(t) model of the present 
5 invention ; 

Fig. 7 is an illustration of an 
operating point; 

Fig. 8 illustrates delay time estimation 
in multiple-stage cells; 
10 Fig. 9 illustrates extraction of circuit 

connection information according to the related 
art; 

Fig. 10 illustrates generation of 
circuit configuration information according to the 
15 related art; 

Fig. 11 illustrates generation of load 
model according to the related art; 

Fig. 12 illustrates logic paths of a 
two-input inverter circuit comprising two 
20 transistors; 

Fig. 13 is a graph showing variation of 
power source voltage with time according to the 
related art; and 

Fig. 14 is a purely capacitive 
25 approximation of a jt-load model subject .to delay 
time estimation. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
First Embodiment 
30 Figs. 1 through 3 show a processing flow 
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for delay time estimation according to a first 
embodiment of the present invention. Fig. 1 shows 
extraction of circuit connection information 
comprising inverters 1 and 2 connected via a wire 3. 
Each of the inverters 1 and 2 is comprised of a 
PMOS transistor and a NMOS transistor (see Fig. 12). 
Circuit configuration information as shown in Fig. 
2 is then generated. A source model 6 corresponding 
to the inverter 1 is modeled by a combination of a 
power source 4 and resistance 5. A load component 8 
comprises a RC distributed constant circuit 7 
corresponding to the wire 3, and input pin 
capacitance Cg corresponding to the inverter 2. 

The input capacitance of the RC 
distributed constant circuit 7 of the load 
component 8 is modeled by capacitance C2. The 
combined output capacitance of the RC distributed 
constant circuit 7 and input pin capacitance Cg is 
modeled by capacitance Cl. The capacitance Cl, 
capacitance C2 and resistance R of the RC 
distributed constant circuit 7 form a rc-load model 
9 as shown in Fig. 3. The description so far 
corresponds to the description of the related art 
with reference to Figs. 9 through 11. 

The present invention introduces a model 
of a power source E , referred to hereinafter as an 
E(t) model, to represent the source model 6. In the 
E(t) model, time-dependence of the source resistor 
of the transistor is represented. As shown in Fig. 
4, a power source voltage waveform is represented 
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using polygonal lines. The voltage varies from 0 to 
V x in a time At x and varies from V L to Vdd in a time 
At 2 . 

The E(t) model according to the 
invention is designed to produce an appropriate 
match with a current flowing from an output pin of 
a gate to a load* Therefore, the current waveform 
is considered in determining the configuration of 
E(t). As shown in Fig. 5, two patterns of current 
waveforms could occur depending on the driving 
capability of the cell and the magnitude of the 
load. In a first pattern, the current makes a 
transition from region 1 in which the current 
increases with time, to region 2 in which the 
current gradually decreases and finally to region 3 
in which the current decreases exponentially. In a 
second pattern, the current makes a transition from 
region 1 in which the current increases with time 
to region 3 in which the current rapidly decreases 
exponentially. 

A consideration is given of the waveform 
of the E(t) model adapted for the current 
characteristic described above. Region 3 is a 
linearity region in which the current determined by 
the time constant of the circuit flows. Therefore, 
E(t)=Vdd is applicable to region 3 . As a 
consequence of this, E(t) functions of different 
waveforms should be applied to regions 1 and 2 
respectively. The E(t) function for region 1 
differs from that of region 2 with respect to 



12 



variation with time- However, both E(t) functions 
for regions 1 and 2 show continuity at points 
corresponding to respective bounds of regions 1 and 
2. As shown in Fig. 4, in the first embodiment, 
straight lines with different gradients are used to 
represent respective functions for respective 
regions to facilitate the ease of computation. 

The time at which E(t) starts to rise 
does not always match a point of time t = 0. The 
delay between the point of time t=0 and the time at 
which E(t) starts to rise is defined as a fixed 
delay t 0 . The waveform of E(t) is expressed as a 
polygonal line showing that the voltage increases 
to El for a first period of time At L after an 
elapse of the fixed delay t 0 . In a subsequent 
period of time At 2 , the voltage increases to Vdd. t 0 
At ir V lf At 2 indicate parameters that describe the 
circuit response. 

The present invention introduces the use 
of a delay library. The delay library is 
constructed by function information At lf V lf At 2 
specifying the Ids-Vds polygonal line at a given 
gate potential and by function information Tslew, 
Tslew being an input slew rate that specifies the 
fixed delay t 0/ and given as a pre-extrac ted delay 
parameter. A large volume of data for the delay 
table according to the related art is not necessary 
The fixed delay t 0 is given as a function of T sl9W . 
Tslew may be defined as a time that elapses before 
the input waveform completes its transition to 0 or 
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Vdd. For example, the fixed delay t 0 is 
approximated by, for example, an equation 
t 0 (Tslew)=K 10 +K 20 *Tslew a . In this case, the function 
information of the input slew rate specifying the 
5 fixed delay and stored in the delay library 

includes K 10 , K 20 and a. 

In pattern 2, region 2 is absent. In 
this case, the waveform of E{t) adapted for 
transition from region 1 to region 3 is given 

10 assuming that At 2 = 0 and V X =E . 

Using the waveform of E(t) shown in Fig. 
4, the gate output v 2 ( t ) of the model at a rising 
edge of the waveform is defined by the following 
equations . 

15 First, under the following definition of 

z, pi, p2 denoting poles of a transfer function, 



RCj RC 2 ^ ^ ^ 

Pi* ^Y^Rq RC2 RSC2 V »n ^4 RRsqcV 



functions f(t,At,V) and a(t,V) are defined as 
20 follows : 

U^^it-^-^^^^-M)!^^"^''^ (2) 



Then, a model waveform description of 
v 2 (t) will now be defined more specifically in the 
following. 
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<Case 1> At x =0 
v 2 (0 



a{ t -tO,V!)+f( t -Iq, At 2 ,E -V t ) 

^ t -tO.Vp+fC t -t Q , At^E -Vp-fC t Atj. Afrj.E -V^ 



<Case 2> At,=0 



v 2 (t)= 



[f( t -to.At^E)-^ t -tQ- 



Atj.A^.E) 



<0s£t£tQ) 

(tQ< tSjtQ+Atj) 

(t^+At^ 

(3) 



(0^ t £tQ) 

(tQ^t^tQ+Atj) 



(4) 



<Case 3> At lf At 3 *0 
v 2 (t)= 



f 0 

f(t -tQ.Atj.Vp 

f( t -tQ.Atj.Vp-fC t -tg- At x . At^Vj) 

+« t -tQ- Atj.Atj.E -Vp 
f( t -Iq. At^Vp-fC t -tQ- At v A^.Vp 



(tQ^t^to+Atj) 
(tQ+At^ t SJtg+At^ At^ 



(t >tg+ Atj+At^ 



+f( t -tQ- Atj,At2»E -Vp-ft t -tQ- Atj- Atj, At2.E - Vj) 



(5) 



A description will now be given of 
determination of the function information At lf . V x 
and At 2 that specify the polygonal lines of Ids-Vds 

10 characteristic. Referring Fig. 5 showing 
transitions of the operating point of the 
transistor, three types of transitions, i.e., 
transition from region 1 to region 2, transition 
from region 1 to region 3 and transition from 

15 region 2 to region 3 are bounded by respective 
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turning points. All of the three turning points 
reside on the Ids-Vds characteristic where Vgs=Vdd. 
This means that the turning points can be defined 
when the Ids-Vds characteristic where Vgs=Vdd is 
available. Region 3 corresponds to the linearity 
region of the MOS transistor, where charging and 
discharging occurs via a fixed resistance Rs , as 
required by the model. Therefore, the Ids-Vds 
characteristic of region 3 should be described 
under the condition Ids=Vds/Rs. Region 2 
corresponds to the saturation region of the MOS 
transistor and described as a region where the 
current gradually decreases. In order to describe 
these features, the Ids-Vds characteristic where 
Vgs=Vdd is represented as shown in Fig. 6, using Rs, 

1 0 and I x . In this model, the region where Vds^RsIi 
is defined as the linearity region, and the region 
where Vds>RsI 1 is defined as the saturation region. 
By defining the saturation current I 0 when Vds=Vdd 
separately, the characteristic where the current 
gradually decreases as Vds decreases (V2 increases) 
in the saturation region is depicted. 

In determining At ir V x and At 2 , Rs , I 0 and 

11 are given as pre-extracted parameters. 

1. Determination of At x 

First, At j, is determined. In determining 
At lf it is assumed that the transistor operates 
according to pattern 2. From the E(t) model circuit 
diagram shown in Fig. 3, it is known that 
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E(t)=v 2 ( t)+Rsi(t ) (6) 

Since we are assuming that pattern 2 takes place, 
5 transition occurs from region 1 to region 3. At a 
turning point t=t 0 +Atl, E(t 0 +Atl)=E so that 

v 2 ( to+At^+Rsif t 0 +At 1 )=E (7) 

10 Since it is assumed that the transistor operates 
according to pattern 2, At 2 =0 so that equation (4) 
is used 

(8) 

Where itt+At^ is a current in the saturation 
15 region. 

(1) According to the simplest transistor model, 
ioc( Vgs-Vth) 2 . Approximating Vgs by a linear 
equation, i<*t 2 . 

(2) If we consider the fixed delay, i=0 at t=t 0 . 
20 (3) Using the linear equation approximation, 

Vgs=Vdd at t=Tslew, resulting in the Ids-Vds 
characteristic of Fig. 6. The current at t=Tslew is 
equal to the current I x at the boundary of the 
linearity region and the saturation region where 
25 Vgs=Vdd. Therefore, 1=^ at t=Tslew. The equation 
that satisfies this condition is given by 
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(9) 



Tslew— t, 



Equation (9) is used for i(t 0 +A 1 ). 

Substituting these equations into (7), we obtain 



1 " i^^? exp( ~ Pl Atl) "^4> exp( ~ P2 A tl)1E +Rsl1 




= 0 



(10) 



5 This equation has a solution at At x . 

When the solution to equation (1) is 



At^Tslew-to, i(t 0 +At 1 )>I 1 so that there is 
digression from the Ids-Vds characteristic of Fig. 
6 since the boundary between region 1 and region 3 

10 must reside on a straight line defining the 

gradient determined by Rs . The reason for this is 
that the initial assumption that pattern 2 is the 
operation pattern failed as a result of transition 
from region 1 to region 2. In this case, we should 

15 proceed assuming that the transistor operates 

according to pattern 1. In pattern 1/ the boundary 
between region 1 and region 2 is at t-to+A^ where 
Vgs=E. Using the same linear equation for Vgs, the 
time is given by t=Tslew so that At^Tslwe-t,. 

20 In an exceptional case, it may be that 

Tslew-t 0 <;0. This is considered as a situation where 
it takes the current to take time before being 
output and where Vgs=E already when the output 
current starts flowing. According to the model, 

25 instant transition to the Ids-Vds where Vgs=Vdd 

occurs. Therefore, a model results in which region 



18 



1 is absent so that the operation starts in region 

2 ( Atj = 0 ) . 

Summarizing the above discussion, At x is 
determined depending on the respective conditions 
5 as described below. 

(1) Case 1: Tslew-t 0 <;0 

At^O 

(2) Case 2: 0<At ^Tslew-t 0 

10 At x satisfies the following equation 



PlP2 Pl<Pl P2> P2 (p I p 2' (Mew-IP) 2 

(11) 



(3) Case 3: when the solution to equation (11) is 

At^Tslew-to 

At^Tslew-to 
15 2. Determination of V x 

The following relationship holds for V x . 

V 2 (t 0 +At x )+Rsi( to+At^Vi (12) 
20 V x is determined based on equation (12). 

(1) Case 1: At x =0 

Vgs=Vdd at t=t 0 so that I(t 0 ) resides on 
the Ids-Vds characteristic of Fig. 6. From the 
25 definition of the fixed delay, v 2 (t 0 )=0 so that 
Vds=Vdd, showing that i(t 0 )=I 0 . Accordingly, 



V^RsIo 



(13) 
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(2) Case : 0<At 1 ^Tslew-t 0 

In this case, the transistor operates 
according to pattern 2 where the transition from 
5 region 1 to region 3 occurs. Therefore, 

V 1= E (14) 

(3) Case 3: At^Tslew-to 

10 Under this condition, the transistor 

operates according to pattern 1 so that Vgs=Vdd at 
t=t 0 +At x . Since 0<v 2 ( t 0 +Ati ) <RsI lf the operating 
point resides at a position shown in Fig- 7. In 
this case, i^o+Ati) is given by the following 

15 equation. 

V2(tn+Att) 

iftQ+At^lQ- Co-Ii) (15) 

Using equation (5), v^to+AtJ is given by 
the following equation. 

(16) 

20 From the equations (12), (15) and (16), 

Vi is given by the following equation. 

RsT 0 

V,= : = : 

1 E— Rsln z Pi~z Pi" 2 1 

1 ~— — {At, -i 7 -exp(-p, At-,) ; -exp(— p? A tj)}-TT 

E-Ralj 1 P1P2 Pl(Pl~P2> P2<Pi~P2> 1 1 At l 

(17) 



3. Determination of At 2 

In case 2, the transistor operates 
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according to pattern 2 so that At 2 is equal to 0. 
In case 1 and case 3, the following equation holds 
at a boundary between region 2 and region 3. 



5 v 2 { 1 0 +At !+At 2 ) +Rs i ( t 0 +At !+At 2 ) =E (18) 

Region 2 is a saturation region where Vgs=Vdd. At 
the boundary between region 2 and region 3, i=l x * 
For v 2 ( t 0 +At 1 +At 2 ) , equation (3) (for case 1) or 
10 equation (for case 3) is used similarly to the 

other processes for determination of the function 
information. 

Therefore, At 2 for case 1 or case 3 

satisfies 

15 

v 2 (t 0 +At 1 +At 2 )+RsI 1 =E ( 19 ) 

In case 1 or case 3, equation (19) has a solution 
at At 2 >0. 

20 A summary is given below how the 

function information is determined for different 
cases . 



<Case 1> Tslew~t 0 i£0 
25 At x =0 

V x =RsI 0 

At 2 should satisfy 
v 2 ( t 0 +At x +At 2 )+RsI x =E 
<Case 2> 0<At 1 ^Tslew-t 0 , where At x is a solution to 
30 the equation below 
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At L 3 

+Rsli - — - ^OSrSI^-fAt] 

{Talew— tQF 

V X =E 

At^O (2 0) 



<Case 3> when the first solution to the equation 

(20) is At a >Tslew-t 0 

A t]l =Tslcw-t 0 

v 

1 E -Rsln z Pl~ z ' P2 -z . _ 1 

v 2 (tgf A ti+ A t2)+Rs Ij = E * W A t2 

(2 1) 

5 The concept described above applies equally to the 
rising edge and the falling edge. 

As described, according to the first 
embodiment, the power source is represeted as two 
types of combinations of straight lines showing 
10 variation with time. Accordingly, the gradual 

decrease of current in the saturation region of the 
transistor is properly reflected so that the delay 
time is estimated in a precise manner. 

15 Second Embodiment 

The method of estimation described above 
applies to a single-stage cell configuration in 
which the input directly controls the gate of the 
output transistor. In a multi-stage cell 

20 configuration (for example, a driving cell) the 
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drain of the transistor having its input controlled 
controls the input of the next transistor so as to 
operate the output transistor. 

The method of estimation for the single- 
5 cell configuration applies equally to the last 

stage of the multi-stage cell. Therefore, for the 
last stage, the same computation as performed for 
the single-cell is performed. In order to effect 
this computation, it is necessary to know the time 

10 to arrive at Vth and the value of Tslew for the 
input of the last stage of the internal node. A 
description will be given with reference to Fig. 8. 

Vth_n denotes the threshold voltage of 
the NMOS transistor in the cell (inverter), Vth_p 

15 denotes the threshold voltage of the PMOS 

transistor in the cell (inverter). TslewP indicates 
the time required for the input waveform of the 
first stage cell to go from Vth_n to the power 
source voltage Vdd. T slewN indicates the time 

20 required for the output waveform of the first stage 
cell (input waveform of the last stage cell) to go 
from Vth_p to the ground voltage 0V. 

The waveform of the first stage may be 
used to determine the time 0 and Tslew of the the 

25 corresponding cell. The reference time for 

computation of the last stage delayed with respect 
to the input to the cell. The delay time is denoted 
t mlt . Tslew of the last stage is different from 
Tslew p . 

30 In the multi-stage cell, the delay time 
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is determined by adding the delay determined using 
Tslew to t mlt . Therefore, in the multi-stage cell, 
the method for estimation of Tslew and T mlt should 
be determined. Both Tslew and t mlt are internal node 
5 values and, therefore, are not considered to be 
dependent on the load connected to the output Y. 
Accordingly, Tslew and t mlt are functions of Tslew p 
specifying the waveform of the input A- The same 
approach employed to estimate the fixed delay of 
10 the single-stage cell is employed to estimate Tslew 
and t mlt . That is, Tslew and t mlt are considered to 
increase as Tslew p increases. Accordingly, the 
internal delay of the multi-stage cell is given by 

Wt^m+W^V* (2 2) 

15 Tslew internal to the multi-stage cell is given by 
Tsfcw=K lt -hK2 t Tstew p a t (2 3) 

where K lm , K 2m , ct m , K lt , K lt and oc t indicate delay 
parameters extracted before the estimation. These 
delay parameters are stored in the delay library. 
20 t nlt and Tslew are determined by using the delay 
library. 

In delay time estimation, computation of 
internal node values is either skipped or not 
skipped depending on whether the single-stage cell 
25 or the multi-stage cell is subject to estimation. 
An efficient approach is to provide the following 
conditions in equations (22) and (23) when the 
single-stage cell is under consideration so that 



24 



t mlt =0 and Tslew=Tslew p . 




(24) 



As described, accordingly to the second embodiment, 
a plurality of logic stages (CMOS logic gates) in a 
5 cell is divided into the last stage and the 

preceding stage(s). By estimating the waveform of 
the gate input at the last stage, the polygonal 
lines of the power source are estimated in a 
precise manner. 

10 

Third Embodiment 



stores a computer program for executing a delay 
time estimation according to the first and second 
15 embodiments. When a computer reads the program from 
the medium and executes the same, the delay time 
estimation of the invention can be efficiently 
performed . 



20 the above-described embodiments, and variations and 
modifications may be made without departing from 
the scope of the present invention. 



A recording medium is provided that 



The present invention is not limited 



to 



